# ENCODING: WINDOWS-1252

library(tidyverse)
library(dplyr)
library(haven)

#Definir ruta
ruta <- "C:/Users/Matias/Dropbox/Proyectos/Fondecyt Participacion Electoral/Articulos/Electoral Winners y Plebiscito/Harvard Dataverse/1. Presidential Election (study 1)/"

# Se carga la base de datos
load(paste0(ruta, "2. Database/MAYO_WIDE.Rdata"))

# PRIMERO SE RECODIFICA LA VARIABLE EDUC DE LA MISMA MANERA EN LA QUE ESTA EN LA BASE DE DATOS DE MAYO. 
# LUEGO SE CALCULA LOS A?OS DE EDUCACION PROMEDIO POR CADA NIVEL EDUCACIONAL EN LA ENCUESTA CASEN, 
# FINALMENTE, DENTRO DE LAS RECODIFICACIONES DE MAYO, SE AGREGA A LA BASE DE DATOS
#casen <- read_dta(paste0(user, "Fondecyt Participacion Electoral/Encuesta Panel Participacion/6_Bases_de_Datos/Casen/Casen 2017.dta"))
#casen <- casen %>% select(esc,educ) %>% filter(educ != 99) %>% 
#  mutate(educ = as.numeric(educ) + 1,
#         rec_educ = case_when(
#           educ == 1 ~ 1,
#           educ == 2 ~ 2,
#           educ == 3 ~ 3,
#           educ == 4 | educ == 5 ~ 4,
#           educ == 6 | educ == 7 ~ 5,
#           educ == 8 ~ 6,
#           educ == 9 ~ 7,
#           educ == 10 ~ 8,
#           educ == 11 | educ == 12 ~ 9,
#           educ == 13 ~ 10)) %>% 
#  group_by(rec_educ) %>% 
#  summarise(prom_anio_educ = round(mean(esc,na.rm = TRUE),2))

# Se cargan datos CASEN para A??os de Educaci??n Promedio
load(file=paste0(ruta,"2. Database/casen_education.Rdata"))

# Recodificacion de variables

mayo_wide <- mayo_wide %>% 
  mutate(
    ## OLA 1
    # VOTO EN OCTUBRE
    part_oct_w01 = factor(P1_w01, 
                          labels = c("No tenia edad para votar", "No vot? en el plebiscito",
                                     "Pens? en votar esta vez, pero no lo hice",
                                     "Normalmente voto, pero esta vez no lo hice",
                                     "Estoy seguro de que vot? en el plebiscito"
                                     )
                          ),
    # VOTO EN OCTUBRE(DICOTOMICA)
    part_octD_w01 = factor(ifelse(
      P1_w01 == 5,1,0),labels = c("No voto en Octubre","Voto en Octubre")
      ),
    # APRUEBO O RECHAZO
    apruebo_w01 = factor(case_when(
      P2_w01 == 1 ~ 2,
      P2_w01 == 2 ~ 0,
      P2_w01 == 3 ~ 1
    ), labels = c("Rechazo","Nulo/Blanco","Apruebo")
    ),
    # APRUEBO O RECHAZO (LOS NULOS Y BLANCOS SON PASADOS A NA)
    aprueboNA_w01 = ifelse(P2_w01 == 3, NA, P2_w01),
    aprueboNA_w01 = factor(ifelse(aprueboNA_w01 == 2,0,1),
                           labels = c("Rechazo","Apruebo")),
    # ELECTORAL WINNER PARA EL PLEBISCITO (NULO/BLANCO COMO DERROTADOS)
    electoral_winner_w01 = factor(ifelse(P2_w01 == 3 | P2_w01 == 1,0,1),
                                  labels = c("Derrotado","Vencedor")),
    # Convenci?n Constitucional o Mixta
    convencion_w01 = factor(case_when(
      P3_w01 == 1 ~ 2,
      P3_w01 == 2 ~ 0,
      P3_w01 == 3 ~ 1
    ), labels = c("Convenci?n Mixta","Nulo/Blanco","Convenci?n Constitucional")
    ),
    # VOTO EN EL 2017
    voto_2017_w01 = factor(P4_w01, 
                           labels = c("No tenia edad para votar",
                                      "No vot? en la elecci?n presidencial",
                                      "Pens? en votar esa vez, pero no lo hice",
                                      "Normalmente voto, pero esa vez no lo hice",
                                      "Estoy seguro de que vote en las elecciones presidenciales"
                                      )
                           ),
    # VOTO EN EL 2017 (DICOTOMICA)
    voto_2017D_w01 = factor(ifelse(
      P4_w01 == 5,1,0), labels = c("No voto en el 2017","Voto en el 2017")
      ),
    # VOTO EN EL 2017 (DISTINGUIENDO A QUIENES NO VOTARON POR EDAD)
    voto_2017T_w01 = factor(case_when(
      P4_w01 == 1 ~ 1,
      P4_w01 >= 2 & P4_w01 <= 4 ~ 2,  
      P4_w01 == 5 ~ 3,
    ), labels = c("Sin Edad Para Votar en 2017","No voto en el 2017","Voto en el 2017")),
    # ELIMINANDO A QUIENES NO TENINAN EDAD PARA VOTAR EN 2017
    voto_2017NA_w01 = factor(ifelse(voto_2017T_w01 == "Sin Edad Para Votar en 2017", NA, voto_2017T_w01),
                             labels = c("No voto en el 2017","Voto en el 2017")),
    # VOTACION EN EL 2017
    candidato_2017_w01 = factor(
      P5_w01, labels = c("Sebastian Pi?era","Alejandro Guillier", "Beatriz S?nchez",
                                           "Jos? Antonio Kast","Marco Enr?quez-Ominami","Carolina Goic",
                                           "Otro", "Nulo/Blanco")
      ),
    # ELECTORAL WINNER - PRESIDENCIAL 2017
    electoral_winner_pdte2017_w01 = factor(ifelse(P5_w01 == 1 | P5_w01 == 4,1,0),
                                            labels = c("Derrotado","Vencedor")), 
    # VOTO EN EL 2016
    voto_2016_w01 = factor(
      P6_w01, labels = c("No tenia edad para votar","No vot? en la elecci?n municipal",
                         "Pens? en votar esa vez pero no lo hice","Normalmente voto, pero esa vez no lo hice",
                         "Estoy seguro de que vot? en las elecciones municipales")),
    # VOTO EN EL 2016 (DICTOMICA)
    voto_2016D_w01 = factor(ifelse(
      P6_w01 == 5,1,0), labels = c("No voto en las municipales","Voto en las municipales")),
    # VOTO EN EL PLEBISCITO DEL 88
    voto_1988_w01 = factor(P7_w01, labels = c("Voto por el S?","Voto por el No",
                                              "No dice el voto","No voto","Sin edad para votar")),
    # QUIENES VOTARON Y NO DECLARAN QU? VOTARON Y QUINES NO VOTARON SON AGRUPADOS EN "OTROS"
    voto_1988rec_w01 = factor(case_when(
      P7_w01 == 1 ~ 0,
      P7_w01 == 2 ~ 1,
      P7_w01 == 3 | P7_w01 == 4 ~ 2,
      P7_w01 == 4 ~ 3)),
    # DISPOSICI?N A IR A VOTAR EL 15 Y 16 DE MAYO
    voto_mayo_w01 = factor(
      P8_w01, labels = c("Definitivamente S?","Probablemente S?",
                         "Probablemente No","Definitivamente No")),
    # DISPOSICI?N A IR A VOTAR EL 15 Y 16 DE MAYO (DICOT?MICA)
    voto_mayoD_w01 = factor(case_when(
      P8_w01 == 1 | P8_w01 == 2 ~ 1,
      P8_w01 == 3 | P8_w01 == 4 ~ 0,
    ), labels = c("No ir? a votar","S? ir? a votar")),
    # POSIBILIDAD DE CONTAGIO COVID
    covid_traslado_w01 = P11_SQ001_w01,
    covid_fila_w01 = P11_SQ002_w01,
    covid_mesa_w01 = P11_SQ003_w01,
    covid_index_w01 = apply(with(mayo_wide, cbind(covid_traslado_w01, covid_fila_w01, covid_mesa_w01)),
                            1, mean, na.rm = TRUE),
    # INDICE DE RIESGO
    riesgo_vida_w01 = P13_SQ001_w01,
    riesgo_financiero_w01 = P13_SQ002_w01,
    riesgo_calle_w01 = P13_SQ003_w01,
    riesgo_salud_w01 = P13_SQ004_w01,
    riesgo_index_w01 = apply(with(mayo_wide,cbind(riesgo_vida_w01, riesgo_financiero_w01,
                                                  riesgo_calle_w01, riesgo_salud_w01)),
                             1, mean, na.rm = TRUE),
    # INTERES EN LA POL?TICA
    pol_interes_w01 = factor(P14_w01, labels = c("Nada interesado","Poco interesado",
                                                  "Algo interesado",
                                                  "Bastante interesado", "Muy interesado")),
    pol_interes3_w01 = factor(case_when(
      P14_w01 <= 2 ~ 0,
      P14_w01 == 3 ~ 1,
      P14_w01 >= 4 ~ 2
    ), labels = c("Nada/Poco interesado","Algo interesado","Bastante/Muy Interesado")),
    # EJE IZQUIERDA-DERECHA
    izquierda_derecha_w01 = factor(case_when(
      P15_SQ001_w01 >= 0 & P15_SQ001_w01 <= 4 ~ 0,
      P15_SQ001_w01 == 5 ~ 1,
      P15_SQ001_w01 >= 6 & P15_SQ001_w01 <= 10 ~ 2,
      P15_SQ001_w01 == 99 ~ 3
    ), labels = c("Izquierda","Centro","Derecha","No Identificado")),
    derecha_w01 = factor(ifelse(izquierda_derecha_w01 == "Derecha",1,0),
                         labels = c("No Derecha","Derecha")),
    izquierda_w01 = factor(ifelse(izquierda_derecha_w01 == "Izquierda",1,0),
                           labels = c("No izquierda","Izquierda")),
    no_id_w01 = factor(ifelse(izquierda_derecha_w01 == "No Identificado",1,0),
                       labels = c("Identificado","No identificado")),
    centro_w01 = factor(ifelse(izquierda_derecha_w01 == "Centro",1,0),
                        labels = c("No centro","Centro")),
    izq_der_cont_w01 = ifelse(P15_SQ001_w01 == 99,NA,P15_SQ001_w01),
    izq_der_cont_w01 = ifelse(is.na(izq_der_cont_w01),
                              round(mean(izq_der_cont_w01,na.rm = TRUE),2),izq_der_cont_w01),
    # CONOCIMIENTO POL?TICO
    pkA_w01 = ifelse(P18A_w01 == 3,1,0),
    pkB_w01 = ifelse(P18B_w01 == 3,1,0),
    pkC_w01 = ifelse(P18C_w01 == 3,1,0),
    pkD_w01 = ifelse(P18D_w01 == 4,1,0),
    pkE_w01 = ifelse(P18E_w01 == 2,1,0),
    pk_w01 = pkA_w01 + pkB_w01 + pkC_w01 + pkD_w01 + pkE_w01,
    # CONFIANZA INTERPERSONAL
    conf1_w01 = P18_SQ001_w01,
    conf2_w01 = P18_SQ002_w01,
    conf3_w01 = P18_SQ003_w01,
    conf_soc_w01 = apply(with(mayo_wide,cbind(conf1_w01,conf2_w01,conf3_w01)),
                         1,mean,na.rm = TRUE),
    # CONFIANZA POLITICA
    trst_gob_w01 = P19_SQ001_w01,
    trst_pp_w01 = P19_SQ002_w01,
    trst_mun_w01 = P19_SQ003_w01,
    trst_just_w01 = P19_SQ004_w01,
    trst_congr_w01 = P19_SQ005_w01,
    trst_servel_w01 = P19_SQ006_w01,
    pol_trst_w01 = apply(with(mayo_wide,cbind(trst_gob_w01,trst_pp_w01,trst_mun_w01,
                                              trst_just_w01,trst_congr_w01,trst_servel_w01)),
                         1,mean,na.rm = TRUE),
    # ACCION POLITICA NO CONVENCIONAL
    marcha_w01 = P20_SQ001_w01,
    caceroleo_w01 = P20_SQ002_w01,
    rrss_w01 = P20_SQ003_w01,
    no_convencional_index_w01 = apply(with(mayo_wide,cbind(marcha_w01,caceroleo_w01,rrss_w01)),
                                      1,mean,na.rm = TRUE),
    # OPINION RESPECTO AL VOTO COMO DEBER
    voto_deber_w01 = P21_SQ001_w01,
    voto_pivot_w01 = P21_SQ002_w01,
    voto_expr_w01 = P21_SQ003_w01,
    voto_supp_w01 = apply(with(mayo_wide,cbind(voto_deber_w01,voto_expr_w01,voto_pivot_w01)),
                          1,mean,na.rm = TRUE),
    # SATISFACCION CON LA DEMOCRACIA
    satis_dem_w01 = factor(P22_w01, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                               "Bastante Satisfecho","Muy Satisfecho")),
    satis_dem3_w01 = factor(case_when(
      P22_w01 <= 2 ~ 0,
      P22_w01 == 3 ~ 1,
      P22_w01 >= 4 ~ 2
    ), labels = c("Insatisfecho","Ni satisfecho ni insatisfecho","Satisfecho")),
    # APROBACION SP
    aprobacion_sp_w01 = factor(P23_w01,labels = c("Aprueba","Desaprueba")),
    # CORRUPCION POLITICA
    corrup_pol_w01 = P24_w01,
    corrup_empr_w01 = P25_w01,
    corrup_index_w01 = apply(with(mayo_wide,cbind(corrup_pol_w01,corrup_empr_w01)),
                             1,mean,na.rm = TRUE),
    # EXPECTATIVAS NUEVA CONSTITUCI?N
    optim_cons1_w01 = P26_SQ001_w01,
    optim_cons2_w01 = 6 - P26_SQ002_w01,
    optim_cons3_w01 = P26_SQ003_w01,
    optim_cons4_w01 = P26_SQ004_w01,
    optim_cons5_w01 = 6 - P26_SQ005_w01,
    optim_cons_index_w01 = apply(with(mayo_wide,cbind(optim_cons1_w01,optim_cons2_w01,optim_cons3_w01,
                                                      optim_cons4_w01,optim_cons5_w01)),
                                 1,mean,na.rm = TRUE),
    # VALORACION ESTALLIDO SOCIAL
    estallido1_w01 = P27_SQ001_w01,
    estallido2_w01 = P27_SQ002_w01,
    estallido3_w01 = 6 - P27_SQ003_w01,
    estallido4_w01 = P27_SQ004_w01,
    estallido5_w01 = 6 - P27_SQ005_w01,
    octubrismo_w01 = apply(with(mayo_wide,cbind(estallido1_w01,
                                                             estallido2_w01,estallido3_w01,
                                                             estallido4_w01,estallido5_w01)),
                                        1,mean,na.rm = TRUE),
    # SOCIOECON?MICOS
    #SEXO
    mujer_w01 = factor(NOMBRE_SEXO_w01),
    #EDAD
    edad_w01 = EDAD_PANELISTA_w01,
    ## NIVEL EDUCACIONAL ENTREVISTADO
    educa_w01 = factor(case_when(
      educ_w01 >= 1 & educ_w01 <= 4 ~ 0,
      educ_w01 == 5 ~ 1,
      educ_w01 == 6 | educ_w01 == 7 ~ 2,
      educ_w01 == 8 ~ 3,
      educ_w01 > 8 & educ_w01 <= 10 ~ 4
    ), labels = c("Secundaria Incomp","Secundaria Comp","Tec. Superior Incp/Comp",
                  "Univer Incomp","Univer Comp")),
    # A?OS DE EDUCACION PROMEDIO PARA CADA NIVEL EDUCACIONAL SEG?N LO HECHO CON LA CASEN
    anio_educ_prom_w01 = case_when(
      educ_w01 == as.numeric(casen[1,1]) ~ as.numeric(casen[1,2]),
      educ_w01 == as.numeric(casen[2,1]) ~ as.numeric(casen[2,2]),
      educ_w01 == as.numeric(casen[3,1]) ~ as.numeric(casen[3,2]),
      educ_w01 == as.numeric(casen[4,1]) ~ as.numeric(casen[4,2]),
      educ_w01 == as.numeric(casen[5,1]) ~ as.numeric(casen[5,2]),
      educ_w01 == as.numeric(casen[6,1]) ~ as.numeric(casen[6,2]),
      educ_w01 == as.numeric(casen[7,1]) ~ as.numeric(casen[7,2]),
      educ_w01 == as.numeric(casen[8,1]) ~ as.numeric(casen[8,2]),
      educ_w01 == as.numeric(casen[9,1]) ~ as.numeric(casen[9,2]),
      educ_w01 == as.numeric(casen[10,1]) ~ as.numeric(casen[10,2])),
    ## SITUACION OCUPACIONAL
    ocup_w01 = factor(ocupa_w01, labels = c("Trabaja de manera remunerada con jornada completa",
                                         "Trabaja de manera remunerada a tiempo parcial o con trabajos ocacionales",
                                         "Estudia y Trabaja", "S?lo estudia","Jubilado","Desempleado",
                                         "Realiza tareas no remuneradas","Inhabilitado para trabajar",
                                         "No estudia, no trabaja, y no busca trabajo")),
    ## TRABAJO ACTUAL
    trab_w01 = factor(trab_ocup_w01, labels = c("Empleado u obrero sector privado","Empleado u obrero sector p?blico",
                                         "Miembro de FFAA y Seguridad","Empleador","Independiente sin empleados",
                                         "Familiar no remunerado","Serivicio Dom?stico")),
    ## TRAMO INGRESO LIQUIDO
    C7_AUX_w01 = ifelse(ingreso_w01 == 99, NA, ingreso_w01),
    ingr_liq_w01 = factor(C7_AUX_w01, labels = c("Menos de $330.000 mensuales",
                                             "$330.001 - $450.000 mensulaes",
                                             "$450.001 - $550.000 mensuales",
                                             "$550.001 - $650.000 mensuales",
                                             "$650.001 - $750.000 mensuales",
                                             "$750.001 - $900.000 mensuales",
                                             "$900.001 - $1.100.000 mensuales",
                                             "$1.100.001 - $1.400.000 mensuales",
                                             "$1.400.001 - $1.800.000 mensuales",
                                             "$1.800.001 - $2.300.000 mensuales",
                                             "$2.300.001 - $2.800.000 mensuales",
                                             "$2.800.001 - $3.800.000 mensuales",
                                             "M?s de 3.800.000 mensuales"
                                             )),
    ## RELIGION
    religion_w01 = factor(ifelse(is.na(relig_w01) & !is.na(relig_other_w01),8,relig_w01),
                          labels = c("Cat?lico","Evang?lico","Protestante","Judio",
                                     "Creyente no Adherente","Ateo","Agn?stico","Otro","Ninguna")),
    prac_relig_w01 = 7 - p_relig_w01,
    ### VER COMO RECODIFICAR PRACTICA RELIGIOSA
    ## OLA 2
    # ASISTENCIA A ELECCION DE MAYO (NO EXISTEN "1. No tenia edad para votar")
    voto_mayo_w02 = factor(P1_B_w02, labels = c("No vot? en la elecci?n del 15 y 16 de Mayo",
                                                "Pens? en votar esta vez, pero no lo hice",
                                                "Normalmente voto, pero esta vez no lo hice",
                                                "Estoy seguro que vot? en las elecciones")),
    voto_mayoD_w02 = factor(case_when(
      P1_B_w02 >= 2 & P1_B_w02 <= 4 ~ 0,
      P1_B_w02 == 5 ~ 1
    ), labels = c("No voto","Voto")),
    # IMPORTANCIA ELECCIONES DE MAYO
    impor_mayo_w02 = factor(P5_B_w02, labels = c("Nada importante","Poco importante",
                                                 "Algo importante","Bastante importante",
                                                 "Muy importante")),
    impor_mayoD_w02 = factor(case_when(
      P5_B_w02 <= 2 ~ 0,
      P5_B_w02 == 3 ~ 1,
      P5_B_w02 >= 4 ~ 2
    ), labels = c("Sin importancia","Algo de imporatncia","Mucha importancia")),
    # INVOLUCRAMIENTO ELECTORAL
    inform_cand_w02 = P6_SQ001_B_w02,
    trab_camp_w02 = P6_SQ002_B_w02,
    rrss_cand_w02 = P6_SQ003_B_w02,
    fam_cand_w02 = P6_SQ004_B_w02,
    firm_cand_w02 = P6_SQ005_B_w02,
    din_cand_w02 = P6_SQ006_B_w02,
    invol_elec_w02 = apply(with(mayo_wide,cbind(inform_cand_w02,trab_camp_w02,rrss_cand_w02,
                                                fam_cand_w02,firm_cand_w02,din_cand_w02)),
                           1,mean,na.rm = TRUE),
    # LISTA CONVENCIONALES CONSTITUYENTES POR LA QUE VOTO
    lista_cc_w02 = factor(P7_B_w02, labels = c("Vamos por Chile","Lista del Apruebo","Apruebo Dignidad",
                                               "Independientes por una Nueva Constituci?n",
                                               "Lista del Pueblo","Otra lista (Partidos Pol?ticos)",
                                               "Otra lista (Independientes)","No recuerda")),
    electoral_winner_cc_w02 = factor(ifelse(P7_B_w02 == 3 | P7_B_w02 == 5,1,0),
                                     labels = c("Derrotado","Vencedor")),
    # PARTIDO POLITICO DEL CANDIDATO A CC VOTADO
    pp_cc_w02 = factor(P8_B_w02, labels = c("PRep","UDI","RN","Evopoli","DC","PR","PPD","PS","PRO",
                                            "PL","Ciudadanos","RD","CS","Comunes","PC","PI","FRVS",
                                            "PH","PEV","No recuerda")),
    electoral_winner_ccpp_w02 = factor(case_when(
      P8_B_w02 <= 7 | P8_B_w02 >= 9 & P8_B_w02 <= 11 | P8_B_w02 >= 16 & P8_B_w02 < 21 ~ 0,
      P8_B_w02 == 8 | P8_B_w02 >= 12 & P8_B_w02 <= 15 ~ 1,
      # QUIENES NO SE ACUERDAN POR QUIEN VOTARON, SON TRATADOS COMO NA
    ), labels = c("Derrotado","Vencedor")),
    # LISTA ALCALDES POR LA QUE VOTO
    lista_alc_w02 = factor(P10_B_w02, labels = c("Republicanos","Chile Vamos","Unidad por el Apruebo",
                                                 "Unidos por la Dignidad","Frente Amplio",
                                                 "Chile Digno, Verde y Soberano","Dignidad Ahora",
                                                 "Otra lista (Partidos Pol?ticos)",
                                                 "Otra lista (Independientes)","No recuerda")),
    # PARTIDO POLITICO CANDIDATO ALCALDE
    pp_alc_w02 = factor(P11_B_w02, labels = c("PRep","UDI","RN","Evopoli","DC","PR","PPD","PS","PRO",
                                              "PL","Ciudadanos","RD","CS","Comunes","PC","PI","FRVS",
                                              "PH","PEV","No recuerda")),
    # LISTA GORE
    lista_gore_w02 = factor(P13_B_w02, labels = c("Republicanos","Chile Vamos","Unidad Constituyente",
                                                  "Frente Amplio","Por Dignidad Regional","Regionalistas Verdes",
                                                  "Humanicemos Chile","Otra lista (Partidos Pol?ticos)",
                                                  "Otra lista (Independientes)","No recuerda")),
    # PARTIDO POL?TICO GORE ("10. PL" no existe en la base de datos)
    pp_gore_w02 = factor(P14_B_w02, labels = c("PRep","UDI","RN","Evopoli","DC","PR","PPD","PS","PRO",
                                               "Ciudadanos","RD","CS","Comunes","PC","PI","FRVS",
                                               "PH","PEV","Candidato Independiente","No recuerda")),
    # RESULTADOS GORE (SOLO QUIENES VOTARON)
    result_gore_w02 = factor(P16_B_w02, labels = c("Sali? electo","Pas? a segunda vuelta","No sali? electo",
                                                   "No sabe")),
    # RESULTADOS ALCALDE (SOLO QUIENES VOTARON)
    result_alc_w02 = factor(P17_B_w02, labels = c("Sali? electo","No sali? electo","No sabe")),
    # RESULTADOS CONVENCIONAL (SOLO QUIENES VOTARON)
    result_cc_w02 = factor(P19_B_w02, labels = c("Sali? electo","No sali? electo", "No sabe")),
    # SATISFACCION RESULTADOS ELECCIONES DE MAYO (SOLO QUIENES VOTARON)
    satis_elec_may_w02 = factor(P25_B_w02, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                                      "Bastante Satisfecho","Muy Satisfecho")),
    satis_elec_may3_w02 = factor(case_when(
      P25_B_w02 <= 2 ~ 0,
      P25_B_w02 == 3 ~ 1,
      P25_B_w02 > 3 & P25_B_w02 <= 5 ~ 2
    ), labels = c("Insatisfecho","Ni Satisfecho ni Insatisfecho","Satisfecho")),
    # RAZONES PARA NO VOTAR VER SI SE PUEDE RECODIFICAR DE ALGUNA MANERA 
    # COMBINACI?N DE LAS PERCEPCIONES DE DESEMPE?O EN UNA SOLA VARIABLE 
    # PARA OBTENER LA PERCPECION DE LAS PERSONAS QUE VOTARON O NO VOTARON, SOLO HAY QUE FILTRAR
    ## CONVENCION CONSTITUCIONAL
    mejor_cc_w02 = factor(coalesce(P22_B_w02,P22b_B_w02), 
                          labels = c("Vamos por Chile","Lista del Apruebo","Apruebo Dignidad",
                                     "Independientes por una Nueva Constituci?n","Lista del Pueblo",
                                     "Otro pacto (Partidos Pol?ticos)", "Otro pacto (Independientes)")),
    ## GOBERNADORES REGIONALES
    mejor_gore_w02 = factor(coalesce(P23_B_w02,P23b_B_w02), 
                            labels = c("Republicanos","Chile Vamos","Unidad Constituyente",
                                       "Frente Amplio","Por Dignidad Regional","Regionalistas Verdes",
                                       "Humanicemos Chile","Independientes")),
    ## ALCALDES
    mejor_alc_w02 = factor(coalesce(P24_B_w02,P24b_B_w02), 
                           labels = c("Republicanos","Chile Vamos","Unidad por el Apruebo",
                                      "Unidad por la Dignidad","Frente Amplio",
                                      "Chile Digno, Verde y Soberano","Dignidad Ahora","Independientes")),
    # INTERES EN POLITICA
    pol_interes_w02 = factor(P14_w02, labels = c("Nada interesado","Poco interesado",
                                                 "Algo interesado",
                                                 "Bastante interesado", "Muy interesado")),
    pol_interes3_w02 = factor(case_when(
      P14_w02 <= 2 ~ 0,
      P14_w02 == 3 ~ 1,
      P14_w02 >= 4 ~ 2
    ), labels = c("Nada/Poco interesado","Algo interesado","Bastante/Muy Interesado")),
    # IDENTIFICACION IZQUIERDA-DERECHA
    izquierda_derecha_w02 = factor(case_when(
      P15_SQ001_w02 >= 0 & P15_SQ001_w02 <= 4 ~ 0,
      P15_SQ001_w02 == 5 ~ 1,
      P15_SQ001_w02 >= 6 & P15_SQ001_w02 <= 10 ~ 2,
      P15_SQ001_w02 == 77 ~ 3
    ), labels = c("Izquierda","Centro","Derecha","No Identificado")),
    derecha_w02 = factor(ifelse(izquierda_derecha_w02 == "Derecha",1,0), 
                         labels = c("No Derecha","Derecha")),
    izquierda_w02 = factor(ifelse(izquierda_derecha_w02 == "Izquierda",1,0), 
                           labels = c("No izquierda","Izquierda")),
    no_id_w02 = factor(ifelse(izquierda_derecha_w02 == "No Identificado",1,0), 
                       labels = c("Identificado","No identificado")),
    centro_w02 = factor(ifelse(izquierda_derecha_w02 == "Centro",1,0), 
                        labels = c("No centro","Centro")),
    izq_der_cont_w02 = ifelse(P15_SQ001_w02 == 77,NA,P15_SQ001_w02),
    izq_der_cont_w02 = ifelse(is.na(izq_der_cont_w02) & atricion_w02 == "Permanece",
                              round(mean(izq_der_cont_w02,na.rm = TRUE),1),izq_der_cont_w02),
    # INDICADORES DE CONFIANZA INTERPERSONAL
    conf1_w02 = P18_SQ001_w02,
    conf2_w02 = P18_SQ002_w02,
    conf3_w02 = P18_SQ003_w02,
    conf_soc_w02 = apply(with(mayo_wide,cbind(conf1_w02,conf2_w02,conf3_w02)),
                         1,mean,na.rm = TRUE),
    # CONFIANZA POLITICA
    trst_gob_w02 = P19_SQ001_w02,
    trst_pp_w02 = P19_SQ002_w02,
    trst_mun_w02 = P19_SQ003_w02,
    trst_just_w02 = P19_SQ004_w02,
    trst_congr_w02 = P19_SQ005_w02,
    trst_servel_w02 = P19_SQ006_w02,
    pol_trst_w02 = apply(with(mayo_wide,cbind(trst_gob_w02,trst_pp_w02,trst_mun_w02,
                                              trst_just_w02,trst_congr_w02,trst_servel_w02)),
                         1,mean,na.rm = TRUE),
    # OPINION RESPECTO AL VOTO
    voto_deber_w02 = P21_SQ001_w02,
    voto_pivot_w02 = P21_SQ002_w02,
    voto_expr_w02 = P21_SQ003_w02,
    voto_supp_w02 = apply(with(mayo_wide,cbind(voto_deber_w02,voto_expr_w02,voto_pivot_w02)),
                          1,mean,na.rm = TRUE),
    # SATISFACCION CON LA DEMOCRACIA
    satis_dem_w02 = factor(P22_w02, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                               "Bastante Satisfecho","Muy Satisfecho")),
    satis_dem3_w02 = factor(case_when(
      P22_w02 <= 2 ~ 0,
      P22_w02 == 3 ~ 1,
      P22_w02 >= 4 ~ 2
    ), labels = c("Insatisfecho","Ni satisfecho ni insatisfecho","Satisfecho")),
    # CORRUPCION POLITICA
    corrup_pol_w02 = P24_w02,
    corrup_empr_w02 = P25_w02,
    corrup_index_w02 = apply(with(mayo_wide,cbind(corrup_pol_w02,corrup_empr_w02)),
                             1,mean,na.rm = TRUE),
    # EXPECTATIVAS NUEVA CONSTITUCI?N
    optim_cons1_w02 = P26_SQ001_w02,
    optim_cons2_w02 = 6 - P26_SQ002_w02,
    optim_cons3_w02 = P26_SQ003_w02,
    optim_cons4_w02 = P26_SQ004_w02,
    optim_cons5_w02 = 6 - P26_SQ005_w02,
    optim_cons_index_w02 = apply(with(mayo_wide,cbind(optim_cons1_w02,optim_cons2_w02,optim_cons3_w02,
                                                      optim_cons4_w02,optim_cons5_w02)),
                                 1,mean,na.rm = TRUE),
    # MATERIALISMO
    materialismo_w02 = apply(with(mayo_wide,cbind(P25b_SQ002_B_w02,P25b_SQ005_B_w02,6-P40_SQ002_B_w02,
                                                  P25b_SQ003_B_w02)),
                             1,mean,na.rm = TRUE),
    # PROGRESISMO
    progresismo_w02 = apply(with(mayo_wide,cbind(P35_SQ004_B_w02,P40_SQ001_B_w02)),1,mean,na.rm = TRUE),
    # VOTO PRESIDENCIAL
    voto_presidencial_w02 = factor(P39_B_w02,labels = c("Jos? Antonio Kast","Joaqu?n Lavin","Mario Desbordes",
                                                        "Sebastian Sichel","Ignacio Briones","Yasna Provoste",
                                                        "Paula Narv?ez","Gabriel Boric","Daniel Jadue",
                                                        "Pamela Jiles","Nulo o Blanco","No votar?a")),
    # SOCIOECONOMICOS
    # SEXO
    mujer_w02 = factor(NOMBRE_SEXO_w02),
    edad_w02 = EDAD_PANELISTA_w02,
    ## JEFE DE HOGAR
    jefe_hogar_w02 = factor(jef_hog_w02, labels = c("S?","No")),
    ## NIVEL EDUCACIONAL PRINCIPAL INGRESO
    educa_jefe_w02 = factor(case_when(
      ed_jefe_w02 >= 1 & ed_jefe_w02 <= 4 ~ 1,
      ed_jefe_w02 == 5 ~ 2,
      ed_jefe_w02 >= 6 & ed_jefe_w02 <= 7 ~ 3,
      ed_jefe_w02 == 8 ~ 4,
      ed_jefe_w02 >= 9 & ed_jefe_w02 <= 10 ~ 5
    ), labels = c("Secundaria Incompleta", "Secundaria Completa", "Tec. Superior Incpleta/Completa",
                  "Universitaria Incompleta", "Universitaria Completa")),
    ## BIENES DEL HOGAR
    internet_w02 = case_when(
      C3_SQ001_w02 == 2 ~ 0,
      C3_SQ001_w02 == 1 ~ 1),
    tv_w02 = case_when(
      C3_SQ002_w02 == 2 ~ 0,
      C3_SQ002_w02 == 1 ~ 1),
    pc_w02 = case_when(
      C3_SQ003_w02 == 2 ~ 0,
      C3_SQ003_w02 == 1 ~ 1),
    secadora_w02 = case_when(
      C3_SQ004_w02 == 2 ~ 0,
      C3_SQ004_w02 == 1 ~ 1),
    auto_w02 = case_when(
      C3_SQ005_w02 == 2 ~ 0,
      C3_SQ005_w02 == 1 ~ 1),
    emp_domes_w02 = case_when(
      C3_SQ006_w02 == 2 ~ 0,
      C3_SQ006_w02 == 1 ~ 1),
    bienes_w02 = apply(with(mayo_wide,cbind(internet_w02,tv_w02,pc_w02,
                                            secadora_w02,auto_w02,emp_domes_w02)),
                   1,mean,na.rm = TRUE),
    ## SISTEMA PREVISIONAL
    prevision_w02 = factor(prev_w02, labels = c("FONASA","CAPREDENA","ISAPRE","Ninguno","No sabe")),
    ## COLEGIO HIJOS O HERMANOS
    cole_hijos_hermanos_w02 = factor(estab_educa_hijo_hmno_w02, labels = c("Liceo","Particular subvencionado","Particular privado",
                                                        "Otro","No sabe")),
    ## ANTECEDENTE COVID
    covid_familia_w02 = factor(covid_fam_w02, labels = c("S?","No","No responde")),
    ### OLA 3
    ## PARTICIPACI?N EN EL PLEBISCITO
    part_oct_w03 = factor(P1_C_w03, labels = c("Vot?","No vot?","Sin edad para votar")),
    ## VOTACI?N EN EL PLEBISCITO
    apruebo_w03 = factor(P2_w03, labels = c("Apruebo","Rechazo","Voto Nulo/Blanco")),
    aprueboNA_w03 = factor(ifelse(P2_w03 == 3,NA,P2_w03),
                           labels = c("Apruebo","Rechazo")),
    ## PARTICIPACI?N EN LAS PRIMARIAS DE JULIO
    prim_julio_w03 = factor(P3_C_w03, labels = c("Vot?","Me abstuve")),
    ## VOTACION EN LAS PRIMARIAS DE JULIO
    cand_prim_jul_w03 = factor(P4_C_w03, labels = c("Gabriel Boric","Daniel Jadue",
                                                    "Sebastian Sichel","Joaqu?n Lavin",
                                                    "Mario Desbordes","Ignacio Briones","Nulo/Blanco")),
    cand_prim_julNA_w03 = factor(ifelse(P4_C_w03 == 99,NA,P4_C_w03),
                                 labels = c("Gabriel Boric","Daniel Jadue",
                                            "Sebastian Sichel","Joaqu?n Lavin",
                                            "Mario Desbordes","Ignacio Briones")),
    ## DERROTADO - VENCEDOR PRIMARIAS
    electoral_winner_prim_w03 = factor(case_when(
      P4_C_w03 == 1 | P4_C_w03 == 3 ~ 1,
      P4_C_w03 == 2 | P4_C_w03 > 3 & P4_C_w03 < 99 ~ 0
    ), labels = c("Derrotado", "Vencedor")),
    ## EN CUAL PRIMARIA PARTICIP?
    part_prim_w03 = factor(case_when(
      P4_C_w03 == 1 | P4_C_w03 == 4 ~ 0,
      P4_C_w03 >= 5 & P4_C_w03 <= 7 ~ 1,
    ), labels = c("Primaria de Izquierda","Primaria de Derecha")),
    ## TIEMPO DEDICADO A PENSAR EN LAS ELECCIONES DE NOVIEMBRE
    pens_electoral_w03 = factor(P5_C_w03, labels = c("Mucho","Bastante","Algo","Poco","Nada")),
    pens_electoral_num_w03 = P5_C_w03,
    ## POSIBILIDADES DE PARTICIPAR EN LAS ELECCIONES DE NOVIEMBRE
    prob_part_nov_w03 = factor(P6_C_w03, labels = c("Definitivamente votar?","Probablemente votar?",
                                                    "Probablemente no votar?","Definitivamente no votar?")),
    prob_part_novD_w03 = factor(ifelse(P6_C_w03 >= 3,0,1),
                                labels = c("No votar?","Votar?")),
    ## DECISI?N PRESIDENCIAL
    decision_pdte_w03 = factor(P7_C_w03, labels = c("S?, s? por cu?l candidato votar?",
                                                         "Tengo algunas alternativas",
                                                         "No s? por cu?l candidato votar?")),
    ## EXPECTATIVAS ELECTORALES POR CANDIDATO
    resultado_boric_w03 = P9_1_C_w03,
    resultado_boric10_w03 = cut(resultado_boric_w03, c(seq(0,50,10),Inf),include.lowest = TRUE),
    resultado_jak_w03 = P9_2_C_w03,
    resultado_jak10_w03 = cut(resultado_jak_w03, c(seq(0,50,10),Inf),include.lowest = TRUE),
    resultado_provoste_w03 = P9_3_C_w03,
    resultado_provoste10_w03 = cut(resultado_provoste_w03, c(seq(0,50,10),Inf),include.lowest = TRUE),
    resultado_sichel_w03 = P9_4_C_w03,
    resultado_sichel10_w03 = cut(resultado_sichel_w03, c(seq(0,50,10),Inf),include.lowst = TRUE),
    resultado_meo_w03 = P9_5_C_w03,
    resultado_meo10_w03 = cut(resultado_meo_w03, c(seq(0,50,10),Inf),include.lowest = TRUE),
    resultado_parisi_w03 = P9_6_C_w03,
    resultado_parisi10_w03 = cut(resultado_parisi_w03, c(seq(0,50,10),Inf),include.lowest = TRUE),
    resultado_artes_w03 = P9_7_C_w03,
    resultado_artes10_w03 = cut(resultado_artes_w03, c(seq(0,50,10),Inf),include.lowest = TRUE),
    ## OPCI?N PRESIDENCIAL
    candidato_2021_w03 = factor(P10_C_w03, labels = c("Gabriel Boric","Jos? Antonio Kast",
                                                      "Yasna Provoste","Sebasti?n Sichel",
                                                      "Marco Enr?quez-Ominami","Franco Parisi",
                                                      "Eduardo Artes","Votar? Nulo","No ir? a votar")),
    ## CANDIDATO POR EL CUAL NO VOTARIA
    no_voto_boric_w03 = factor(case_when(
      P11_1_C_w03 == 1 ~ 0,
      P11_1_C_w03 == 0 ~ 1,
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_jak_w03 = factor(case_when(
      P11_2_C_w03 == 1 ~ 0,
      P11_2_C_w03 == 0 ~ 1
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_provoste_w03 = factor(case_when(
      P11_3_C_w03 == 1 ~ 0,
      P11_3_C_w03 == 0 ~ 1
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_sichel_w03 = factor(case_when(
      P11_4_C_w03 == 1 ~ 0,
      P11_4_C_w03 == 0 ~ 1
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_meo_w03 = factor(case_when(
      P11_5_C_w03 == 1 ~ 0,
      P11_5_C_w03 == 0 ~ 1
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_parisi_w03 = factor(case_when(
      P11_6_C_w03 == 1 ~ 0,
      P11_6_C_w03 == 0 ~ 1
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_artes_w03 = factor(case_when(
      P11_7_C_w03 == 1 ~ 0,
      P11_7_C_w03 == 0 ~ 1
    ), labels = c("No votar?a por","Votar?a por")),
    no_voto_nr_w03 = factor(case_when(
      P11_99_C_w03 == 0 ~ 0,
      P11_99_C_w03 == 1 ~ 1,
      is.na(P11_99_C_w03) & atricion_w03 != "Atrici?n" ~ 1
    ), labels = c("No mencionado","Mencionado")),
    ## NUMERO DE CANDIDATOS POR LOS CUALES NO VOTAR?A
    n_no_alternativa_w03 = ifelse(atricion_w03 != "Atrici?n",apply(.[699:705],1,function(x)sum(x == 1)),
                            NA),
    n_alternativa_w03 = 7 - n_no_alternativa_w03,
    covid_traslado_w03 = P11_SQ001_w03,
    covid_fila_w03 = P11_SQ002_w03,
    covid_mesa_w03 = P11_SQ003_w03,
    covid_index_w03 = apply(with(mayo_wide, cbind(covid_traslado_w03, covid_fila_w03, covid_mesa_w03)),
                            1, mean, na.rm = TRUE),
    ## ?QU? ES EL VOTO?
    votar_w03 = factor(P13A_C_w03, labels = c("Un deber","Una opci?n","No estoy seguro")),
    voto_deber_w03 = factor(P13B_C_w03, labels = c("Muy frecuentemente","Algo frecuentemente",
                                                   "No muy frecuentemente")),
    pol_interes_w03 = factor(P14_w03, labels = c("Nada interesado","Poco interesado","Algo interesado",
                                                 "Bastante interesado","Muy interesado")),
    pol_interes3_w03 = factor(case_when(
      P14_w03 == 1 | P14_w03 == 2 ~ 0,
      P14_w03 == 3 ~ 1,
      P14_w03 == 4 | P14_w03 == 5 ~ 2), labels = c("Nada/Poco interesado","Algo interesado",
                                                   "Bastante/Muy Interesado")),
    # EJE IZQUIERDA-DERECHA
    izquierda_derecha_w03 = factor(case_when(
      P15_SQ001_w03 >= 0 & P15_SQ001_w03 <= 4 ~ 0,
      P15_SQ001_w03 == 5 ~ 1,
      P15_SQ001_w03 >= 6 & P15_SQ001_w03 <= 10 ~ 2,
      P15_SQ001_w03 == 99 ~ 3
    ), labels = c("Izquierda","Centro","Derecha","No Identificado")),
    derecha_w03 = factor(ifelse(izquierda_derecha_w03 == "Derecha",1,0),
                         labels = c("No Derecha","Derecha")),
    izquierda_w03 = factor(ifelse(izquierda_derecha_w03 == "Izquierda",1,0),
                           labels = c("No izquierda","Izquierda")),
    no_id_w03 = factor(ifelse(izquierda_derecha_w03 == "No Identificado",1,0),
                       labels = c("Identificado","No identificado")),
    centro_w03 = factor(ifelse(izquierda_derecha_w03 == "Centro",1,0),
                        labels = c("No centro","Centro")),
    izq_der_cont_w03 = ifelse(P15_SQ001_w03 == 99,NA,P15_SQ001_w03),
    izq_der_cont_w03 = ifelse(is.na(izq_der_cont_w03) & atricion_w03 == "Permanece",
                              round(mean(izq_der_cont_w03,na.rm = TRUE),2),izq_der_cont_w03),
    # EJE IZQUIERDA-DERECHA PARA LOS CANDIDATOS
    id_kast_w03 = factor(case_when(
      P15_1_C_w03 >= 0 & P15_1_C_w03 <= 4 ~ 0,
      P15_1_C_w03 == 5 ~ 1,
      P15_1_C_w03 >= 6 & P15_1_C_w03 <= 10 ~ 2,
      P15_1_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_kast_cont_w03 = ifelse(P15_1_C_w03 == 99,NA,P15_1_C_w03),
    id_kast_cont_w03 = ifelse(is.na(id_kast_cont_w03),
                              round(mean(id_kast_cont_w03,na.rm = TRUE),2),id_kast_cont_w03),
    no_sabe_kast_w03 = factor(ifelse(P15_1_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    id_provoste_w03 = factor(case_when(
      P15_2_C_w03 >= 0 & P15_2_C_w03 <= 4 ~ 0,
      P15_2_C_w03 == 5 ~ 1,
      P15_2_C_w03 >= 6 & P15_2_C_w03 <= 10 ~ 2,
      P15_2_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_provoste_cont_w03 = ifelse(P15_2_C_w03 == 99,NA,P15_2_C_w03),
    id_provoste_cont_w03 = ifelse(is.na(id_provoste_cont_w03),
                              round(mean(id_provoste_cont_w03,na.rm = TRUE),2),id_provoste_cont_w03),
    no_sabe_provoste_w03 = factor(ifelse(P15_2_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    id_boric_w03 = factor(case_when(
      P15_3_C_w03 >= 0 & P15_3_C_w03 <= 4 ~ 0,
      P15_3_C_w03 == 5 ~ 1,
      P15_3_C_w03 >= 6 & P15_3_C_w03 <= 10 ~ 2,
      P15_3_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_boric_cont_w03 = ifelse(P15_3_C_w03 == 99,NA,P15_3_C_w03),
    id_boric_cont_w03 = ifelse(is.na(id_boric_cont_w03),
                              round(mean(id_boric_cont_w03,na.rm = TRUE),2),id_boric_cont_w03),
    no_sabe_boric_w03 = factor(ifelse(P15_3_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    id_sichel_w03 = factor(case_when(
      P15_4_C_w03 >= 0 & P15_4_C_w03 <= 4 ~ 0,
      P15_4_C_w03 == 5 ~ 1,
      P15_4_C_w03 >= 6 & P15_4_C_w03 <= 10 ~ 2,
      P15_4_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_sichel_cont_w03 = ifelse(P15_4_C_w03 == 99,NA,P15_4_C_w03),
    id_sichel_cont_w03 = ifelse(is.na(id_sichel_cont_w03),
                              round(mean(id_sichel_cont_w03,na.rm = TRUE),2),id_sichel_cont_w03),
    no_sabe_sichel_w03 = factor(ifelse(P15_4_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    id_meo_w03 = factor(case_when(
      P15_5_C_w03 >= 0 & P15_5_C_w03 <= 4 ~ 0,
      P15_5_C_w03 == 5 ~ 1,
      P15_5_C_w03 >= 6 & P15_5_C_w03 <= 10 ~ 2,
      P15_5_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_meo_cont_w03 = ifelse(P15_5_C_w03 == 99,NA,P15_5_C_w03),
    id_meo_cont_w03 = ifelse(is.na(id_meo_cont_w03),
                              round(mean(id_meo_cont_w03,na.rm = TRUE),2),id_meo_cont_w03),
    no_sabe_meo_w03 = factor(ifelse(P15_5_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    id_artes_w03 = factor(case_when(
      P15_6_C_w03 >= 0 & P15_6_C_w03 <= 4 ~ 0,
      P15_6_C_w03 == 5 ~ 1,
      P15_6_C_w03 >= 6 & P15_6_C_w03 <= 10 ~ 2,
      P15_6_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_artes_cont_w03 = ifelse(P15_6_C_w03 == 99,NA,P15_6_C_w03),
    id_artes_cont_w03 = ifelse(is.na(id_artes_cont_w03),
                              round(mean(id_artes_cont_w03,na.rm = TRUE),2),id_artes_cont_w03),
    no_sabe_artes_w03 = factor(ifelse(P15_6_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    id_parisi_w03 = factor(case_when(
      P15_7_C_w03 >= 0 & P15_7_C_w03 <= 4 ~ 0,
      P15_7_C_w03 == 5 ~ 1,
      P15_7_C_w03 >= 6 & P15_7_C_w03 <= 10 ~ 2,
      P15_7_C_w03 > 10 ~ 3
    ), labels = c("Izquieda","Centro","Derecha","No Sabe")),
    id_parisi_cont_w03 = ifelse(P15_7_C_w03 == 99,NA,P15_7_C_w03),
    id_parisi_cont_w03 = ifelse(is.na(id_parisi_cont_w03),
                              round(mean(id_parisi_cont_w03,na.rm = TRUE),2),id_parisi_cont_w03),
    no_sabe_parisi_w03 = factor(ifelse(P15_7_C_w03 == 99,1,0),labels = c("No Identifica","Identifica")),
    # CONOCIMIENTO POL?TICO
    pkA_w03 = ifelse(P17A_C_w03 == 2,1,0),
    pkB_w03 = ifelse(P17B_C_w03 == 2,1,0),
    pkC_w03 = ifelse(P18B_w03 == 3,1,0),
    pkD_w03 = ifelse(P17D_C_w03 == 3,1,0),
    pkE_w03 = ifelse(P17E_C_w03 == 4,1,0),
    pk_w03 = pkA_w03 + pkB_w03 + pkC_w03 + pkD_w03 + pkE_w03,
    # APROBACION 4to RETIRO 10%
    apruebo_4retiro_w03 = factor(P18_C_w03, labels = c("Muy en contra","En contra",
                                                       "Ni a favor ni en contra","A favor",
                                                       "Muy a favor")),
    apruebo_4retiro_rec_w03 = factor(case_when(
      P18_C_w03 == 1 | P18_C_w03 == 2 ~ 0,
      P18_C_w03 == 3 ~ 1,
      P18_C_w03 == 4 | P18_C_w03 == 5 ~ 2
    ), labels = c("En contra","Ni a favor ni en contra","A favor")),
    # CIVISMO
    civismo_1_w03 = P19_1_C_w03,
    civismo_2_w03 = P19_2_C_w03,
    civismo_3_w03 = P19_3_C_w03,
    civismo_4_w03 = P19_4_C_w03,
    civismo_index_w03 = apply(with(mayo_wide, cbind(civismo_1_w03, civismo_2_w03,
                                                    civismo_3_w03, civismo_4_w03)),
                              1, mean, na.rm = TRUE),
    # CONFIANZA INTERPERSONAL
    conf1_w03 = P18_SQ001_w03,
    conf2_w03 = P18_SQ002_w03,
    conf3_w03 = P18_SQ003_w03,
    conf_soc_w03 = apply(with(mayo_wide,cbind(conf1_w01,conf2_w01,conf3_w01)),
                         1,mean,na.rm = TRUE),
    # CONFIANZA INSTITUCIONAL
    trst_gob_w03 = P19_SQ001_w03,
    trst_pp_w03 = P19_SQ002_w03,
    trst_mun_w03 = P19_SQ003_w03,
    trst_just_w03 = P19_SQ004_w03,
    trst_congr_w03 = P19_SQ005_w03,
    trst_servel_w03 = P19_SQ006_w03,
    trst_conv_w03 = P19_SQ007_w03,
    ## (SIN CONVENCION)
    pol_trst_w03 = apply(with(mayo_wide,cbind(trst_gob_w03,trst_pp_w03,trst_mun_w03,
                                              trst_just_w03,trst_congr_w03,trst_servel_w03)),
                         1,mean,na.rm = TRUE),
    ## (CON CONVENCION)
    pol_trst2_w03 = apply(with(mayo_wide,cbind(trst_gob_w03,trst_pp_w03,trst_mun_w03,
                                       trst_just_w03,trst_congr_w03,trst_servel_w03,
                                       trst_conv_w03)),
                          1,mean,na.rm = TRUE),
    # OPINION VOTO
    voto_deber_w03 = P21_SQ001_w03,
    voto_pivot_w03 = P21_SQ002_w03,
    voto_expr_w03 = P21_SQ003_w03,
    voto_supp_w03 = apply(with(mayo_wide,cbind(voto_deber_w03,voto_expr_w03,voto_pivot_w03)),
                          1,mean,na.rm = TRUE),
    # SATISFACCION CON LA DEMOCRACIA
    satis_dem_w03 = factor(P22_w03, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                               "Bastante Satisfecho","Muy Satisfecho")),
    satis_dem3_w03 = factor(case_when(
      P22_w03 <= 2 ~ 0,
      P22_w03 == 3 ~ 1,
      P22_w03 >= 4 ~ 2
    ), labels = c("Insatisfecho","Ni satisfecho ni insatisfecho","Satisfecho")),
    # CORRUPCION EN LA POLITICA
    corrup_pol_w03 = P24_w03,
    corrup_empr_w03 = P25_w03,
    corrup_index_w03 = apply(with(mayo_wide,cbind(corrup_pol_w03,corrup_empr_w03)),
                             1,mean,na.rm = TRUE),
    # OPTIMISMO CONSTITUCIONAL
    optim_cons1_w03 = P26_SQ001_w03,
    optim_cons2_w03 = 6 - P26_SQ002_w03,
    optim_cons3_w03 = P26_SQ003_w03,
    optim_cons4_w03 = P26_SQ004_w03,
    optim_cons5_w03 = 6 - P26_SQ005_w03,
    optim_cons_index_w03 = apply(with(mayo_wide,cbind(optim_cons1_w03,optim_cons2_w03,optim_cons3_w03,
                                                      optim_cons4_w03,optim_cons5_w03)),
                                 1,mean,na.rm = TRUE),
    # OCTUBRISMO
    estallido1_w03 = P27_SQ001_w03,
    estallido2_w03 = P27_SQ002_w03,
    estallido3_w03 = 6 - P27_SQ003_w03,
    estallido4_w03 = P27_SQ004_w03,
    estallido5_w03 = 6 - P27_SQ005_w03,
    octubrismo_w03 = apply(with(mayo_wide,cbind(estallido1_w03,
                                                             estallido2_w03,estallido3_w03,
                                                             estallido4_w03,estallido5_w03)),
                                        1,mean,na.rm = TRUE),
    ## SOCIODEMOGRAFICOS
    mujer_w03 = factor(sexo_w03, labels = c("Hombre","Mujer")),
    edad_w03 = EDAD_PANELISTA_w03,
    educa_w03 = factor(case_when(
      educ_w03 >= 1 & educ_w03 <= 4 ~ 0,
      educ_w03 == 5 ~ 1,
      educ_w03 == 6 | educ_w03 == 7 ~ 2,
      educ_w03 == 8 ~ 3,
      educ_w03 > 8 & educ_w03 <= 10 ~ 4
    ), labels = c("Secundaria Incomp","Secundaria Comp","Tec. Superior Incp/Comp",
                  "Univer Incomp","Univer Comp")),
    # A?OS DE EDUCACION PROMEDIO PARA CADA NIVEL EDUCACIONAL SEG?N LO HECHO CON LA CASEN
    anio_educ_prom_w03 = case_when(
      educ_w03 == as.numeric(casen[1,1]) ~ as.numeric(casen[1,2]),
      educ_w03 == as.numeric(casen[2,1]) ~ as.numeric(casen[2,2]),
      educ_w03 == as.numeric(casen[3,1]) ~ as.numeric(casen[3,2]),
      educ_w03 == as.numeric(casen[4,1]) ~ as.numeric(casen[4,2]),
      educ_w03 == as.numeric(casen[5,1]) ~ as.numeric(casen[5,2]),
      educ_w03 == as.numeric(casen[6,1]) ~ as.numeric(casen[6,2]),
      educ_w03 == as.numeric(casen[7,1]) ~ as.numeric(casen[7,2]),
      educ_w03 == as.numeric(casen[8,1]) ~ as.numeric(casen[8,2]),
      educ_w03 == as.numeric(casen[9,1]) ~ as.numeric(casen[9,2]),
      educ_w03 == as.numeric(casen[10,1]) ~ as.numeric(casen[10,2])),
    ## OCUPACION
    ocup_w03 = factor(ocupa_w03, labels = c("Trabaja de manera remunerada con jornada completa",
                                            "Trabaja de manera remunerada a tiempo parcial o con trabajos ocacionales",
                                            "Estudia y Trabaja", "S?lo estudia","Jubilado","Desempleado",
                                            "Realiza tareas no remuneradas","Inhabilitado para trabajar",
                                            "No estudia, no trabaja, y no busca trabajo")),
    ## TRABAJO ACTUAL O PASADO
    trab_w03 = factor(trab_ocup_w01, labels = c("Empleado u obrero sector privado","Empleado u obrero sector p?blico",
                                              "Miembro de FFAA y Seguridad","Empleador","Independiente sin empleados",
                                              "Familiar no remunerado","Serivicio Dom?stico")),
    ## TRAMO INGRESO LIQUIDO
    C7_AUX_w03 = ifelse(ingreso_w03 == 99, NA, ingreso_w03),
    ingr_liq_w03 = factor(C7_AUX_w03, labels = c("Menos de $330.000 mensuales", #1
                                                 "$330.001 - $450.000 mensulaes", #2
                                                 "$450.001 - $550.000 mensuales", #3
                                                 "$550.001 - $650.000 mensuales", #4
                                                 "$650.001 - $750.000 mensuales", #5
                                                 "$750.001 - $900.000 mensuales", #6
                                                 "$900.001 - $1.100.000 mensuales", #7
                                                 "$1.100.001 - $1.400.000 mensuales", #8
                                                 "$1.400.001 - $1.800.000 mensuales", #9
                                                 "$1.800.001 - $2.300.000 mensuales", #10
                                                 "$2.300.001 - $2.800.000 mensuales", #11
                                                 "$2.800.001 - $3.800.000 mensuales", #12
                                                 "M?s de 3.800.000 mensuales")), #13
    # NUMERO DE PERSONAS EN EL HOGAR
    n_personas_hog_w03 = factor(ifelse(n_personas_hogar_w03 > 6,7,n_personas_hogar_w03),
                                labels = c("1","2","3","4",
                                           "5","6","7 o m?s")),
    ingr_n_hog_w03 = case_when(
      # N?mero de personas en el hogar == 1
      C7_AUX_w03 == 1 & n_personas_hogar_w03 == 1 ~ 3,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 == 1 ~ 4,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 == 1 ~ 4,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 == 1 ~ 4,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 == 1 ~ 5,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 == 1 ~ 5,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 == 1 ~ 5,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 == 1 ~ 6,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 == 1 ~ 6,
      C7_AUX_w03 >= 10 & n_personas_hogar_w03 == 1 ~ 7,
      # N?mero de personas en el hogar == 2
      C7_AUX_w03 == 1 & n_personas_hogar_w03 == 2 ~ 2,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 == 2 ~ 3,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 == 2 ~ 3,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 == 2 ~ 4,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 == 2 ~ 4,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 == 2 ~ 4,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 == 2 ~ 4,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 == 2 ~ 5,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 == 2 ~ 5,
      C7_AUX_w03 == 10 & n_personas_hogar_w03 == 2 ~ 6,
      C7_AUX_w03 == 11 & n_personas_hogar_w03 == 2 ~ 6,
      C7_AUX_w03 >= 12 & n_personas_hogar_w03 == 2 ~ 7,
      # N?mero de personas en el hogar == 3
      C7_AUX_w03 == 1 & n_personas_hogar_w03 == 3 ~ 2,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 == 3 ~ 2,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 == 3 ~ 3,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 == 3 ~ 3,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 == 3 ~ 3,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 == 3 ~ 4,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 == 3 ~ 4,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 == 3 ~ 4,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 == 3 ~ 5,
      C7_AUX_w03 == 10 & n_personas_hogar_w03 == 3 ~ 5,
      C7_AUX_w03 == 11 & n_personas_hogar_w03 == 3 ~ 6,
      C7_AUX_w03 == 12 & n_personas_hogar_w03 == 3 ~ 6,
      C7_AUX_w03 == 13 & n_personas_hogar_w03 == 3 ~ 7,
      # N?mero de personas en el hogar == 4
      C7_AUX_w03 == 1 & n_personas_hogar_w03 == 4 ~ 1,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 == 4 ~ 2,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 == 4 ~ 2,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 == 4 ~ 3,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 == 4 ~ 3,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 == 4 ~ 3,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 == 4 ~ 4,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 == 4 ~ 4,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 == 4 ~ 4,
      C7_AUX_w03 == 10 & n_personas_hogar_w03 == 4 ~ 5,
      C7_AUX_w03 == 11 & n_personas_hogar_w03 == 4 ~ 5,
      C7_AUX_w03 == 12 & n_personas_hogar_w03 == 4 ~ 6,
      C7_AUX_w03 == 13 & n_personas_hogar_w03 == 4 ~ 7,
      # N?mero de personas en el hogar == 5
      C7_AUX_w03 == 1 & n_personas_hogar_w03 == 5 ~ 1,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 == 5 ~ 2,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 == 5 ~ 2,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 == 5 ~ 2,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 == 5 ~ 3,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 == 5 ~ 3,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 == 5 ~ 3,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 == 5 ~ 4,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 == 5 ~ 4,
      C7_AUX_w03 == 10 & n_personas_hogar_w03 == 5 ~ 5,
      C7_AUX_w03 == 11 & n_personas_hogar_w03 == 5 ~ 5,
      C7_AUX_w03 == 12 & n_personas_hogar_w03 == 5 ~ 6,
      C7_AUX_w03 == 13 & n_personas_hogar_w03 == 5 ~ 7,
      # N?mero de personas en el hogar == 6
      C7_AUX_w03 == 1 & n_personas_hogar_w03 == 6 ~ 1,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 == 6 ~ 1,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 == 6 ~ 2,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 == 6 ~ 2,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 == 6 ~ 2,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 == 6 ~ 3,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 == 6 ~ 3,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 == 6 ~ 4,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 == 6 ~ 4,
      C7_AUX_w03 == 10 & n_personas_hogar_w03 == 6 ~ 4,
      C7_AUX_w03 == 11 & n_personas_hogar_w03 == 6 ~ 5,
      C7_AUX_w03 == 12 & n_personas_hogar_w03 == 6 ~ 6,
      C7_AUX_w03 == 13 & n_personas_hogar_w03 == 6 ~ 7,
      # N?mero de personas en el hogar >= 7
      C7_AUX_w03 == 1 & n_personas_hogar_w03 >= 7 ~ 1,
      C7_AUX_w03 == 2 & n_personas_hogar_w03 >= 7 ~ 1,
      C7_AUX_w03 == 3 & n_personas_hogar_w03 >= 7 ~ 2,
      C7_AUX_w03 == 4 & n_personas_hogar_w03 >= 7 ~ 2,
      C7_AUX_w03 == 5 & n_personas_hogar_w03 >= 7 ~ 2,
      C7_AUX_w03 == 6 & n_personas_hogar_w03 >= 7 ~ 2,
      C7_AUX_w03 == 7 & n_personas_hogar_w03 >= 7 ~ 3,
      C7_AUX_w03 == 8 & n_personas_hogar_w03 >= 7 ~ 3,
      C7_AUX_w03 == 9 & n_personas_hogar_w03 >= 7 ~ 4,
      C7_AUX_w03 == 10 & n_personas_hogar_w03 >= 7 ~ 4,
      C7_AUX_w03 == 11 & n_personas_hogar_w03 >= 7 ~ 5,
      C7_AUX_w03 == 12 & n_personas_hogar_w03 >= 7 ~ 5,
      C7_AUX_w03 == 13 & n_personas_hogar_w03 >= 7 ~ 6),
    ### OLA 4
    #PARTICIPACI?N EN LAS ELECCIONES DE NOVIEMBRE (TODOS LOS PARTICIPANTES TENIAN EDAD PARA VOTAR)
    part_nov_w04 = factor(P1_D_w04, labels = c("No vot? en la elecci?n del 21 de Noviembre",
                                              "Pens? en votar esta vez, pero no lo hice","Normalmente voto, pero esta vez no lo hice",
                                              "Estoy seguro de que vot? en la elecci?n")),
    part_novD_w04 = factor(ifelse(P1_D_w04 < 5,0,1), labels = c("No voto", "Vot?")),
    #VOTACI?N EN EL MISMO LOCAL QUE MAYO DE 2021
    mismo_local_mayo_w04 = factor(case_when(
      P3A_D_w04 == 2 | P3A_D_w04 == 3 ~ 0,
      P3A_D_w04 == 1 ~ 1,
      P3A_D_w04 == 4 ~ 2
    ), labels = c("No vot? en el mismo local","Vot? en el mismo local",
                  "No particip? de las elecciones de Mayo")),
    #VOTACION EN EL MISMO LOCAL QUE MAYO (QUIENES NO VOTARON EN MAYO COMO NA)
    mismo_local_mayoNA_w04 = factor(case_when(
      P3A_D_w04 == 2 | P3A_D_w04 == 3 ~ 0,
      P3A_D_w04 == 1 ~ 1
    ), labels = c("No vot? en el mismo local","Vot? en el mismo local")),
    #DISTANCIA DEL NUEVO LOCAL DE VOTACION
    dist_nuevo_local_w04 = factor(P4a_B_w04, labels = c("M?s cerca","M?s lejos","M?s o menos igual")),
    #LOCAL M?S LEJOS
    local_lejos_w04 = factor(ifelse(P4a_B_w04 == 2,1,0)),
    #LOCAL M?S CERCA
    local_cerca_w04 = factor(ifelse(P4a_B_w04 == 1,1,0)),
    #LOCAL MISMO DISTANCIA
    local_misma_dist_w04 = factor(ifelse(P4a_B_w04 == 3,1,0)),
    # IMPORTANCIA ELECCIONES DE NOVIEMBRE
    impor_nov_w04 = factor(P5_B_w04, labels = c("Nada importante","Poco importante",
                                                 "Algo importante","Bastante importante",
                                                 "Muy importante")),
    impor_mayoD_w04 = factor(case_when(
      P5_B_w04 <= 2 ~ 0,
      P5_B_w04 == 3 ~ 1,
      P5_B_w04 >= 4 ~ 2
    ), labels = c("Sin importancia","Algo de imporatncia","Mucha importancia")),
    # INVOLUCRAMIENTO ELECTORAL NOVIEMBRE
    inform_cand_w04 = P6_SQ001_B_w04,
    trab_camp_w04 = P6_SQ002_B_w04,
    rrss_cand_w04 = P6_SQ003_B_w04,
    fam_cand_w04 = P6_SQ004_B_w04,
    firm_cand_w04 = P6_SQ005_B_w04,
    din_cand_w04 = P6_SQ006_B_w04,
    invol_elec_w04 = apply(with(mayo_wide,cbind(inform_cand_w02,trab_camp_w02,rrss_cand_w02,
                                                fam_cand_w02,firm_cand_w02,din_cand_w02)),
                           1,mean,na.rm = TRUE),
    # VOTO PRESIDENCIAL
    voto_pres_w04 = factor(P9_D_w04, labels = c("Gabriel Boric","Jos? Antonio Kast",
                                               "Yasna Provoste","Sebasti?n Sichel",
                                               "Eduardo Artes","Marco Enriquez-Ominami",
                                               "Franco Parisi","No me acuerdo")),
    # VOTO PRESIDENCIAL (NO ME ACUERDO COMO NA)
    voto_presNA_w04 = factor(ifelse(voto_pres_w04 == "No me acuerdo",NA,voto_pres_w04),
                             labels = c("Gabriel Boric","Jos? Antonio Kast",
                                        "Yasna Provoste","Sebasti?n Sichel",
                                        "Eduardo Artes","Marco Enriquez-Ominami",
                                        "Franco Parisi")),
    # VENCEDOR - DERROTADOS (NO ME ACUERDO COMO NA)
    electoral_winner_prim_vuel2021_w04 = factor(case_when(
      P9_D_w04 == 1 | P9_D_w04 == 2 ~ 1,
      P9_D_w04 > 2 & P9_D_w04 < 99 ~ 0,
    )), 
    # VOTO PRESIDENCIAL SEG?N IZQUIERDA-DERECHA (PARISI COMO SIN ID; NO ME ACUERDO COMO NA)
    voto_presID_w04 = factor(case_when(
      P9_D_w04 == 1 | P9_D_w04 == 3 | P9_D_w04 == 5 | P9_D_w04 == 5 ~ 1,
      P9_D_w04 == 2 | P9_D_w04 == 4 ~ 2,
      P9_D_w04 == 6 ~ 3
    ), labels = c("Izquierda","Derecha","Sin Id")),
    # VOTO LISTA DIPUTADOS
    list_dip_w04 = factor(ifelse(P10_D_w04 == 99,NA,P10_D_w04), 
                          labels = c("Chile Podemos M?s","Partido de la Gente",
                                     "Nuevo Pacto Social","Uni?n Patri?tica",
                                     "Dignidad Ahora","Frente Social Cristiano",
                                     "Apruebo Dignidad","Otra lista de partidos",
                                     "Independientes")),
    # VOTO LISTA DIPUTADOS SEG?N IZQUIERDA-DERECHA (PARTIDO DE LA GENTE, OTRA LISTA E INDEPENDIENTES COMO SIN ID.)
    list_dipID_w04 = factor(case_when(
      P10_D_w04 == 3 | P10_D_w04 == 4 | P10_D_w04 == 5 | P10_D_w04 == 7 ~ 1,
      P10_D_w04 == 1 | P10_D_w04 == 6 ~ 2,
      P10_D_w04 == 2 | P10_D_w04 >= 8 ~ 3
    ), labels = c("Izquierda","Derecha","Sin Id")),
    # VOTO PARTIDOS DIPUTADO
    part_dip_w04 = factor(ifelse(P11_D_w04 == 99,21,P11_D_w04),
                          labels = c("Partido Rep?blicano","Uni?n Dem?crata Independiente",
                                     "Renovaci?n Nacional","Evoluci?n Pol?tica","Democracia Cristiana",
                                     "Partido Radical","Partido por la Democracia","Partido Socialista",
                                     "Partido Progresista","Partido Liberal","Ciudadanos",
                                     "Revoluci?n Democr?tica","Convergencia Social","Comunes",
                                     "Partido Comunista", "Partido Igualdad",
                                     "Federaci?n Regionalista Verde Social","Partido Humanista",
                                     "Partido Ecologista Verde","Partido de la Gente","No me acuerdo")),
    # VENCEDOR - DERROTADO ELECCIONES PARLAMENTARIAS (NO ME ACUERDO SON DEJADOS COMO NA)
    electoral_winner_dip_w04 = factor(case_when(
      P11_D_w04 >= 1 & P11_D_w04 <= 3 ~ 1,
      P11_D_w04 == 8 | P11_D_w04 >= 12 & P11_D_w04 <= 15 ~ 1,
      P11_D_w04 == 20 ~ 1,
      P11_D_w04 >= 4 & P11_D_w04 <= 7 ~ 0,
      P11_D_w04 >= 9 & P11_D_w04 <= 11 ~ 0,
      P11_D_w04 > 15 & P11_D_w04 < 20 ~ 0,
    ), labels = c("Derrotado", "Vencedor")),
    # VOTO PARTIDOS DIPUTADO (NO ME ACUERDO COMO NA)
    part_dipNA_w04 = factor(ifelse(P11_D_w04 == 99,NA,P11_D_w04), labels = c("Partido Rep?blicano","Uni?n Dem?crata Independiente",
                                                                               "Renovaci?n Nacional","Evoluci?n Pol?tica","Democracia Cristiana",
                                                                               "Partido Radical","Partido por la Democracia","Partido Socialista",
                                                                               "Partido Progresista","Partido Liberal","Ciudadanos",
                                                                               "Revoluci?n Democr?tica","Convergencia Social","Comunes","Partido Comunista",
                                                                               "Partido Igualdad","Federaci?n Regionalista Verde Social",
                                                                               "Partido Humanista","Partido Ecologista Verde","Partido de la Gente")),
    # OPCI?N PRESIDENCIAL: ELECTO O NO 
    elect_pres_w04 = factor(P16_D_w04, labels = c("Sali? electo","Pas? a segunda vuelta","No sali? electo",
                                                    "No sabe")),
    # OPCI?N DIPUTADO: ELECTO O NO
    elect_dip_w04 = factor(P17_D_w04, labels = c("Sali? electo","No salio electo","No sabe")),
    # LISTA OBTUVO REPRESENTACI?N A NIVEL NACIONAL: SI O NO
    list_repr_w04 = factor(P18_D_w04, labels = c("S?, obtuvo al menos un esca?o","No obtuvo esca?o","No sabe")),
    # PACTO CON MEJOR APOYO ELECTORAL
    mejor_dipVOT_w04 = factor(P19_D_w04, labels = c("Chile Podemos M?s","Partido de la Gente","Nuevo Pacto Social",
                                                   "Frente Social Cristiano","Apruebo Dignidad")),
    # SATISFACCI?N CON RESULTADO ELECTORAL
    satis_elec_w04 = factor(P20_D_w04, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                                      "Bastante Satisfecho","Muy Satisfecho")),
    satis_elec_w04 = factor(case_when(
      P20_D_w04 <= 2 ~ 0,
      P20_D_w04 == 3 ~ 1,
      P20_D_w04 > 3 & P20_D_w04 <= 5 ~ 2
    ), labels = c("Insatisfecho","Ni Satisfecho ni Insatisfecho",
                  "Satisfecho")),
    # COMBINACI?N DE LAS PERCEPCIONES DE DESEMPE?O EN UNA SOLA VARIABLE 
    # PARA OBTENER LA PERCPECION DE LAS PERSONAS QUE VOTARON O NO VOTARON, SOLO HAY QUE FILTRAR
    mejor_dip_w04 = factor(coalesce(P19_D_w04,P22_D_w04), 
                            labels = c("Chile Podemos M?s","Partido de la Gente","Nuevo Pacto Social",
                                       "Frente Social Cristiano","Apruebo Dignidad")),
    # INTERES EN POL?TICA
    pol_interes_w04 = factor(P14_w04, labels = c("Nada interesado","Poco interesado",
                                                 "Algo interesado",
                                                 "Bastante interesado", "Muy interesado")),
    pol_interes3_w04 = factor(case_when(
      P14_w04 <= 2 ~ 0,
      P14_w04 == 3 ~ 1,
      P14_w04 >= 4 ~ 2
    ), labels = c("Nada/Poco interesado","Algo interesado","Bastante/Muy Interesado")),
    # EJE IZQUIERDA-DERECHA
    izquierda_derecha_w04 = factor(case_when(
      P15_SQ001_w04 >= 0 & P15_SQ001_w04 <= 4 ~ 0,
      P15_SQ001_w04 == 5 ~ 1,
      P15_SQ001_w04 >= 6 & P15_SQ001_w04 <= 10 ~ 2,
      P15_SQ001_w04 == 99 ~ 3
    ), labels = c("Izquierda","Centro","Derecha","No Identificado")),
    derecha_w04 = factor(ifelse(izquierda_derecha_w01 == "Derecha",1,0),
                           labels = c("No Derecha","Derecha")),
    izquierda_w04 = factor(ifelse(izquierda_derecha_w01 == "Izquierda",1,0),
                           labels = c("No izquierda","Izquierda")),
    no_id_w04 = factor(ifelse(izquierda_derecha_w01 == "No Identificado",1,0),
                       labels = c("Identificado","No identificado")),
    centro_w04 = factor(ifelse(izquierda_derecha_w01 == "Centro",1,0),
                        labels = c("No centro","Centro")),
    izq_der_cont_w04 = ifelse(P15_SQ001_w04 == 99,NA,P15_SQ001_w04),
    izq_der_cont_w04 = ifelse(is.na(izq_der_cont_w01) & atricion_w04 == "Permanece",
                            round(mean(izq_der_cont_w04,na.rm = TRUE),2),izq_der_cont_w04),
    # SIMPATIA PARTIDARIA
    simp_part_w04 = factor(P16_w04,labels = c("S?","No")),
    partido_simpatizado_w04 = factor(P17_w04, labels = c("Partido Republicano","Uni?n Dem?crata Independiente",
                                                           "Renovaci?n Nacional","Evoluci?n Pol?tica","Democracia Cristiana",
                                                           "Partido Radical","Partido por la Democracia","Partido Socialista",
                                                           "Revoluci?n Democr?tica","Convergencia Social","Comunes",
                                                           "Partido Comunista","Partido Humanista","Partido de la Gente")),
    grado_simpatia_w04 = factor(P30_B_w04, labels = c("Poco cercano","Algo cercano",
                                                        "Bastante cercano","Muy cercano")),
    # CONFIANZA INTERPERSONAL
    conf1_w04 = P18_SQ001_w04,
    conf2_w04 = P18_SQ002_w04,
    conf3_w04 = P18_SQ003_w04,
    conf_soc_w04 = apply(with(mayo_wide,cbind(conf1_w04,conf2_w04,conf3_w04)),
                         1,mean,na.rm = TRUE),
    # CONFIANZA POLITICA
    trst_gob_w04 = P19_SQ001_w04,
    trst_pp_w04 = P19_SQ002_w04,
    trst_mun_w04 = P19_SQ003_w04,
    trst_just_w04 = P19_SQ004_w04,
    trst_congr_w04 = P19_SQ005_w04,
    trst_servel_w04 = P19_SQ006_w04,
    trst_conv_w04 = P19_SQ007_w04,
    ## (SIN CONVENCION)
    pol_trst_w03 = apply(with(mayo_wide,cbind(trst_gob_w03,trst_pp_w03,trst_mun_w03,
                                              trst_just_w03,trst_congr_w03,trst_servel_w03)),
                         1,mean,na.rm = TRUE),
    ## (CON CONVENCION)
    pol_trst2_w03 = apply(with(mayo_wide,cbind(trst_gob_w03,trst_pp_w03,trst_mun_w03,
                                               trst_just_w03,trst_congr_w03,trst_servel_w03,
                                               trst_conv_w03)),
                          1,mean,na.rm = TRUE),
    # OPINION RESPECTO AL VOTO
    voto_deber_w04 = P21_SQ001_w04,
    voto_pivot_w04 = P21_SQ002_w04,
    voto_expr_w04 = P21_SQ003_w04,
    voto_supp_w04 = apply(with(mayo_wide,cbind(voto_deber_w04,voto_expr_w04,
                                               voto_pivot_w04)),
                          1,mean,na.rm = TRUE),
    # SATISFACCI?N CON LA DEMOCRACIA
    satis_dem_w04 = factor(P22_w04, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                                 "Bastante Satisfecho","Muy Satisfecho")),
    satis_dem3_w04 = factor(case_when(
      P22_w04 <= 2 ~ 0,
      P22_w04 == 3 ~ 1,
      P22_w04 >= 4 ~ 2
    ), labels = c("Insatisfecho","Ni satisfecho ni insatisfecho","Satisfecho")),
    # CORRUPCION EN LA POLITICA
    corrup_pol_w04 = P24_w04,
    corrup_empr_w04 = P25_w04,
    corrup_index_w04 = apply(with(mayo_wide,cbind(corrup_pol_w04,corrup_empr_w04)),
                             1,mean,na.rm = TRUE),
    # OPTIMISMO CONSTITUCIONAL
    optim_cons1_w04 = P26_SQ001_w04,
    optim_cons2_w04 = 6 - P26_SQ002_w04,
    optim_cons3_w04 = P26_SQ003_w04,
    optim_cons4_w04 = P26_SQ004_w04,
    optim_cons5_w04 = 6 - P26_SQ005_w04,
    optim_cons_index_w04 = apply(with(mayo_wide,cbind(optim_cons1_w04,optim_cons2_w04,optim_cons3_w04,
                                                      optim_cons4_w04,optim_cons5_w04)),
                                 1,mean,na.rm = TRUE),
    # VOTACION SEGUNDA VUELTA
    votacion_segunda_w04 = factor(P39_D_w04, labels =c("Jos? Antonio Kast","Gabriel Boric",
                                                       "No lo he decidido a?n","Votar?a nulo o blanco",
                                                       "No votar?a")),
    ## SOCIODEMOGR?FICOS
    ## SEXO
    mujer_w04 = factor(CODIGO_SEXO_w04),
    edad_w04 = EDAD_PANELISTA_w04,
    jefe_hogar_w04 = factor(C1_w04, labels = c("S?","No")),
    ## NIVEL EDUCACIONAL PRINCIPAL INGRESO
    educa_jefe_w04 = factor(case_when(
      ed_jefe_w04 >= 1 & ed_jefe_w04 <= 4 ~ 1,
      ed_jefe_w04 == 5 ~ 2,
      ed_jefe_w04 >= 6 & ed_jefe_w04 <= 7 ~ 3,
      ed_jefe_w04 == 8 ~ 4,
      ed_jefe_w04 >= 9 & ed_jefe_w04 <= 10 ~ 5
    ), labels = c("Secundaria Incompleta", "Secundaria Completa", "Tec. Superior Incpleta/Completa",
                  "Universitaria Incompleta", "Universitaria Completa")),
    ## BIENES DEL HOGAR
    internet_w04 = case_when(
      C3_SQ001_w04 == 2 ~ 0,
      C3_SQ001_w04 == 1 ~ 1),
    tv_w04 = case_when(
      C3_SQ002_w04 == 2 ~ 0,
      C3_SQ002_w04 == 1 ~ 1),
    pc_w04 = case_when(
      C3_SQ003_w04 == 2 ~ 0,
      C3_SQ003_w04 == 1 ~ 1),
    secadora_w04 = case_when(
      C3_SQ004_w04 == 2 ~ 0,
      C3_SQ004_w04 == 1 ~ 1),
    auto_w04 = case_when(
      C3_SQ005_w04 == 2 ~ 0,
      C3_SQ005_w04 == 1 ~ 1),
    emp_domes_w04 = case_when(
      C3_SQ006_w04 == 2 ~ 0,
      C3_SQ006_w04 == 1 ~ 1),
    bienes_w04 = apply(with(mayo_wide,cbind(internet_w04,tv_w04,pc_w04,
                                            secadora_w04,auto_w04,emp_domes_w04)),
                   1,mean,na.rm = TRUE),
    ## SISTEMA PREVISIONAL
    prevision_w04 = factor(prev_w04, labels = c("FONASA","CAPREDENA","ISAPRE","Ninguno","No sabe")),
    ## ANTECEDENTE COVID
    covid_familia_w04 = factor(covid_fam_w04, labels = c("S?","No","No responde")),
    ### OLA 5 (SEGUNDA VUELTA PRESIDENCIAL)
    part_1v_w05 = factor(P1_E_w05, labels = c("No vot? en la elecci?n del 19 de diciembre", "Pens? en votar esta vez, pero no lo hice",
                                              "Normalmente voto, pero esta vez no lo hice", "Estoy seguro que vot? en la elecci?n del 19 de diciembre")),
    part_1vD_w05 = factor(ifelse(P1_E_w05 < 5,0,1), labels = c("No vot?", "Vot?")),
    # IMPORTANCIA ELECCIONES DE LAS ELECCIONES PRESIDENCIALES
    impor_presi_w05 = factor(P2_E_w05, labels = c("Nada importante","Poco importante",
                                                 "Algo importante","Bastante importante",
                                                 "Muy importante")),
    impor_presiD_w05 = factor(case_when(
      P2_E_w05 <= 2 ~ 0,
      P2_E_w05 == 3 ~ 1,
      P2_E_w05 >= 4 ~ 2
    ), labels = c("Sin importancia","Algo de imporatncia","Mucha importancia")),
    ## INVOLUCRAMIENTO ELECTORAL
    inform_cand_w05 = P6_SQ001_B_w05,
    trab_camp_w05 = P6_SQ002_B_w05,
    rrss_cand_w05 = P6_SQ003_B_w05,
    fam_cand_w05 = P6_SQ004_B_w05,
    invol_elec_w05 = apply(with(mayo_wide,cbind(inform_cand_w05,trab_camp_w05,rrss_cand_w05,
                                                fam_cand_w05)),
                           1,mean,na.rm = TRUE),
    ## CONVERSACION SOBRE ELECCIONES
    miemb_fam_w05 = P4_1_E_w05,
    amigos_w05 = P4_2_E_w05,
    vecinos_w05 = P4_3_E_w05,
    trabajo_w05 = P4_4_E_w05,
    prom_conv_w05 = apply(with(mayo_wide,cbind(miemb_fam_w05, amigos_w05,
                                               vecinos_w05, trabajo_w05)),
                          1, mean, na.rm = TRUE),
    # APOYOS DE GRUPO
    ap_miemb_fam_w05 = factor(P5_1_E_w05, labels = c("Otro Candidato","Indeciso","Apoyaban al mismo candidato")),
    ap_amigos_w05 = factor(P5_2_E_w05, labels = c("Otro cansadidato","Indeciso","Apoyaban al mismo candidato")),
    ap_vecinos_w05 = factor(P5_3_E_w05, labels = c("Otro candidato","Indeciso","Apoyaban al mismo candidato")),
    ap_trabajo_w05 = factor(P5_4_E_w05, labels = c("Otro candidato","Indeciso","Apoyaban al mismo candidato")),
    ## VOTO EN LA SEGUNDA VUELTA
    voto_2vuelta_w05 = factor(P7_E_w05, labels = c("Gabriel Boric","Jos? Antonio Kast",NA)),
    ## VENCEDOR - DERROTADO ELECTORAL
    electoral_winner_pdte2021_w05 = factor(P7_E_w05, levels = c(2,1,3),
                                           labels = c("Derrotado","Vencedor",NA)),
    ## SATISFACCION RESULTADOS EN LAS ELECCIONES
    satis_2v_w05 = factor(P8_E_w05, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                                      "Bastante Satisfecho","Muy Satisfecho")),
    satis_2v_w05 = factor(case_when(
      P8_E_w05 <= 2 ~ 0,
      P8_E_w05 == 3 ~ 1,
      P8_E_w05 > 3 & P8_E_w05 <= 5 ~ 2
    ), labels = c("Insatisfecho","Ni Satisfecho ni Insatisfecho","Satisfecho")),
    ## SENTIMIENTO SEGUNDA VUELTA
    contento_w05 = P9_1_E_w05,
    preocupado_w05 = P9_2_E_w05,
    orgulloso_w05 = P9_3_E_w05,
    triste_w05 = P9_4_E_w05,
    enojado_w05 = P9_5_E_w05,
    esperanzado_w05 = P9_6_E_w05,
    indiferente_w05 = P9_7_E_w05,
    sent_positivo_w05 = contento_w05 + (6 - preocupado_w05) + orgulloso_w05 + 
      (6-triste_w05) + (6-enojado_w05) + esperanzado_w05,
    ## VER RECODIFICAR RAZONES NO VOTAR
    ## INTERES EN POLITICA
    pol_interes_w05 = factor(P14_w05, labels = c("Nada interesado","Poco interesado",
                                                 "Algo interesado",
                                                 "Bastante interesado", "Muy interesado")),
    pol_interes3_w05 = factor(case_when(
      P14_w05 <= 2 ~ 0,
      P14_w05 == 3 ~ 1,
      P14_w05 >= 4 ~ 2
    ), labels = c("Nada/Poco interesado","Algo interesado","Bastante/Muy Interesado")),
    # EJE IZQUIERDA-DERECHA
    izquierda_derecha_w05 = factor(case_when(
      P15_SQ001_w05 >= 0 & P15_SQ001_w05 <= 4 ~ 0,
      P15_SQ001_w05 == 5 ~ 1,
      P15_SQ001_w05 >= 6 & P15_SQ001_w05 <= 10 ~ 2,
      P15_SQ001_w05 == 99 ~ 3
    ), labels = c("Izquierda","Centro","Derecha","No Identificado")),
    derecha_w05 = factor(ifelse(izquierda_derecha_w05 == "Derecha",1,0),
                         labels = c("No Derecha","Derecha")),
    izquierda_w05 = factor(ifelse(izquierda_derecha_w05 == "Izquierda",1,0),
                           labels = c("No izquierda","Izquierda")),
    no_id_w05 = factor(ifelse(izquierda_derecha_w05 == "No Identificado",1,0),
                       labels = c("Identificado","No identificado")),
    centro_w05 = factor(ifelse(izquierda_derecha_w05 == "Centro",1,0),
                        labels = c("No centro","Centro")),
    izq_der_cont_w05 = ifelse(P15_SQ001_w05 == 99,NA,P15_SQ001_w05),
    izq_der_cont_w05 = ifelse(is.na(izq_der_cont_w05) & atricion_w05 == "Permanece",
                              round(mean(izq_der_cont_w05,na.rm = TRUE),2),izq_der_cont_w05),
    # SIMPATIA PARTIDARIA
    simp_part_w05 = factor(P16_w05,labels = c("S?","No")),
    partido_simpatizado_w05 = factor(P17_w05, labels = c("Partido Republicano","Uni?n Dem?crata Independiente",
                                                         "Renovaci?n Nacional","Evoluci?n Pol?tica","Democracia Cristiana",
                                                         "Partido Radical","Partido por la Democracia","Partido Socialista",
                                                         "Revoluci?n Democr?tica","Convergencia Social","Comunes",
                                                         "Partido Comunista","Partido Humanista","Partido de la Gente")),
    grado_simpatia_w05 = factor(P30_B_w05, labels = c("Poco cercano","Algo cercano",
                                                      "Bastante cercano","Muy cercano")),
    # CONFIANZA INTERPERSONAL
    conf1_w05 = P18_SQ001_w05,
    conf2_w05 = P18_SQ002_w05,
    conf3_w05 = P18_SQ003_w05,
    conf_soc_w05 = apply(with(mayo_wide,cbind(conf1_w05,conf2_w05,conf3_w05)),
                         1,mean,na.rm = TRUE),
    # CONFIANZA POLITICA
    trst_gob_w05 = P19_SQ001_w05,
    trst_pp_w05 = P19_SQ002_w05,
    trst_just_w05 = P19_SQ004_w05,
    trst_congr_w05 = P19_SQ005_w05,
    trst_servel_w05 = P19_SQ006_w05,
    pol_trst_w05 = apply(with(mayo_wide,cbind(trst_gob_w05,trst_pp_w05,
                                              trst_just_w05,trst_congr_w05,trst_servel_w05)),
                         1,mean,na.rm = TRUE),
    # OPINION RESPECTO AL VOTO COMO DEBER
    voto_deber_w05 = P21_SQ001_w05,
    voto_pivot_w05 = P21_SQ002_w05,
    voto_expr_w05 = P21_SQ003_w05,
    voto_supp_w05 = apply(with(mayo_wide,cbind(voto_deber_w05,voto_expr_w05,
                                               voto_pivot_w05)),
                          1,mean,na.rm = TRUE),
    # SATISFACCION CON LA DEMOCRACIA
    satis_dem_w05 = factor(P22_w05, labels = c("Nada Satisfecho","Poco Satisfecho","Algo Satisfecho",
                                               "Bastante Satisfecho","Muy Satisfecho")),
    satis_dem3_w05 = factor(case_when(
      P22_w05 <= 2 ~ 0,
      P22_w05 == 3 ~ 1,
      P22_w05 >= 4 ~ 2
    ), labels = c("Insatisfecho","Ni satisfecho ni insatisfecho","Satisfecho")),
    # CORRUPCION POLITICA
    corrup_pol_w05 = P24_w05,
    corrup_empr_w05 = P25_w05,
    corrup_index_w05 = apply(with(mayo_wide,cbind(corrup_pol_w05,corrup_empr_w05)),
                             1,mean,na.rm = TRUE),
    mujer_w05 = factor(CODIGO_SEXO_w05),
    edad_w05 = EDAD_PANELISTA_w05,
    # IDENTIFICADOR DE ENCUESTA
    id_encuesta_w01 = 1,
    id_encuesta_w02 = ifelse(atricion_w02 == "Permanece",1,0),
    id_encuesta_w03 = ifelse(atricion_w03 == "Permanece",1,0),
    id_encuesta_w04 = ifelse(atricion_w04 == "Permanece",1,0),
    id_encuesta_w05 = ifelse(atricion_w05 == "Permanece",1,0),
    id_encuesta = paste0(id_encuesta_w01,id_encuesta_w02,id_encuesta_w03,
                         id_encuesta_w04,id_encuesta_w05,sep = ""))

## SE GUARDA LA BASE DE DATOS.
save(mayo_wide, file = paste0(ruta, "2. Database/MAYO_WIDE_RECODE.Rdata"))
